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Quantization of Thermodynamics, Supersecondary 
Quantization, and a New Variational Principle 
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oo . V. P. Maslov* 
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In solving the problem of finding a temperature distribution which, at zero temper- 
ature, corresponds to superfluidity, i.e., to nonzero energy, the author tried to quantize 
free energy. This was done on the basis of supersecondary quantization whose spe- 
rial case is the usual secondary quantization for bosons and with the help of which 
new representations of the Schrodinger equation were obtained. The supersecondary 
quantization allowed the author to construct a variational method whose zero approx- 
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Abstract 



1.1 Different pictures of N bodies and their representations 



imation are the Hartree-Fock and Bogolyubov-BCSch variational principles. This 
method works especially well in the case of not a large number of particles. The new 
quantization and the variational method are of general character and can be used in 
the quantum field theory. 

> 

oo ■ 

o> : 1 Introduction 

<N 

o 

There are two aspects in the many-body problem. The many bodies can be located in our 
three-dimensional space and move in some way. They can collide or interact. All this takes 
place in our three-dimensional space. However, we can consider these bodies in a different 
way. Namely, we can treat these A bodies as one point (here the term "A bodies" means "A 
material points" ) in the three-dimensional space. In this three-dimensional space this single 
point moves in some way, and we can observe its motion. In some cases one of these two 
approaches (two aspects) is efficient. In other cases the other approach can be useful. So, 
from the mathematical viewpoint, it is simpler to deal with a single point. From the physical 
viewpoint, since we live in the three-dimensional space, we would like to understand how 
material points behave in our three-dimensional space. This is a more natural "picture." 
For example, in the three-dimensional space we can study a particle located in the field of 
all other particles. In particular, this approach underlies Vlasov's theory of self-consistent 
field. In this theory it is assumed that all particles are approximately the same and have the 
same distribution. We can assume that one particle behaves as if it interacts with another 
particle, and this other particle is the same particle, i.e., it has the same distribution as the 
first particle. 
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It turns out that this approach can be realized not only approximately, as it was realized 
by Vlasov, but exactly by using the secondary quantization, i.e., in other words, in the 
three-dimensional space we consider the behavior of operators rather than the behavior of 
functions. Namely, we study some creation and annihilation operators. For the first time, 
this was realized by Dirac in quantum mechanics. However, it is not necessary to consider 
just the quantum mechanics. So, Schoenberg studied the classical mechanics (approximately 
in 1953-1956) and applied this secondary quantization method to classical objects and to 
classical statistical physics. 

Indeed, this is a very general thing. For example, I together with my pupil used this 
method in the case of N fields. We acted as if it was the third quantization, which, in 
fact, is a version of the secondary quantization [1]. Here we have a natural generalization. 
If we consider iV particles as a single particle in the three-dimensional space, why cannot 
we consider, for example, any two particles of these N particles? I can treat this pair 
as one particle in the six-dimensional space and then study these pairs of six-dimensional 
particles. Each of them is similar to one particle in the 3iV/2-dimensional space. So I can 
simultaneously consider a pair of particles and a single particle. In other words, there exist 
some six-dimensional particles (pairs) and some three-dimensional particles and all this is 
embedded in the corresponding space. If there are k pairs and m "single" particles, then we 
have 2k + m = 3N, where N is the general number of particles. 

In fact, in addition to the coordinate of the particle (the location of the point), there 
can be some other degrees of freedom. For example, in the classical mechanics there can be 
momenta. In the quantum theory, there can be spins and some other degrees of freedom. 

In the theory developed here, I introduce one more degree of freedom, the minimal degree 
which is just a number. In addition to the location of a particle in some part of the space, this 
particle is equipped with a number. This viewpoint had a good interpretation in the early 
quantum mechanics. In the Bohr theory we deal with electrons. We can interchange their 
places, but nothing changes. This is the identity principle. However, besides these electrons, 
there are orbit numbers. They can also be treated as numbers assigned to electrons. 

For example, I have N particles and I consider some pairs of particles and single particles, 
which I have numbered. Perhaps, it is convenient also to number the pairs. Whether I have 
numbered then or not, this does not change anything, since this numbering is not related 
to the fact that I fix these particles and they cannot interchange their places. Thus all this 
is the same picture. These are simply different representations of the same phenomenon. 
The secondary quantization allows us to write equations in the three-dimensional space 
but with operators, i.e., an operator equation. This operator equation is equivalent to one 
equation in the 3n-dimensional space (which is simply an equation) not for operators but 
for functions. The same problem is obtained in terms of representations, which I propose 
and call supersecondary quantized. 

1.2 Series and serials 

Probabilities of transition from level to level play an essential role. In particular, as is known, 
the energy quantums issued by an atom correspond to transitions from level to level. If the 
probability of transition is zero or very small, i.e., if this transition is practically forbidden, 
then the corresponding levels are treated as if they belong to different series. In fact, the 



2 



partition into series means a partition into systems that are isolated from each other. If 
a given series of energy levels contains the minimal level, which is the principal state of 
the series, then this is a metastable state. If the system is in this state and transitions 
to the lower levels of other series are almost forbidden, then the system stays in this state 
sufficiently long. 

The existence of a superfluid liquid or a superconducting current corresponds to such 
metastable series. The current can flow without friction during 100000 years, which means 
that the series is almost stable. In a finite large volume, where we consider a system with a 
large number of particles, the velocities and currents are discrete, and each value determines 
its own metastable series. A set of series that differ in velocities we call a serial. 

To each picture discussed above, there corresponds its own serial in the asymptotic ap- 
proximation. Although the physicists do not say this directly, in reality, they try to find a 
serial such that the lower level of its lower series (i.e., the level at the zero velocity) would 
coincide with the principal state of the entire original system. 

In the above-mentioned paper by N. N. Bogolyubov (1947), such a serial was found. It 
corresponds to the following picture: N particles are in the three-dimensional space and all 
their numbers coincide. 

However, it turns out and this will be proved in the present work that the serial cor- 
responding to the picture in which pairs of particles are not numbered but the particles 
are numbered has a level which is lower than that of the Bogolyubov serial. In this sense, 
apparently, the serial corresponds to the true phenomenon of superfluidity and, to a larger 
extent, corresponds to the experiment. Although, if the nature "set" a system to another 
serial, the system will stay in this state for a long time. 

However, it is possible that resonances may occur between different series. In this case the 
series is destroyed. An example is provided by the one-dimensional Schrodinger equation 
with two minima. If the maximum between the minima is very large and if we have a 
semiclassical approximation, then to each hollow there corresponds its own series and if 
there are no resonances, then transitions /from one hollow to the other are almost forbidden 
(the "tunnel" effects). If some eigenvalues of two series coincide, then the resonance is 
possible and in this case the position of a particle in one of the hollows at this eigenvalue 
is destroyed. For example, in the potential symmetric with respect to the maximum, all 
eigenvalues "almost" coincide, all are in resonance, and the system does not split into series. 

1.3 Quantization of free energy 

We assume that the eigenfunctions of the Schrodinger equation (the Hamilton operator) 
remain the same. From each eigenvalue we subtract the temperature multiplied by entropy 
corresponding to this eigenvalue. In the simplest case of a discrete spectrum, as is known [2] , 
this is the logarithm of the multiplicity of the eigenvalue. 

As is known, the Hamilton operator can be represented as the sum of its eigenvalues 
multiplied by the projectors on the corresponding subspaces of eigenfunctions. The free 
energy operator, as was already pointed out, is the sum of the same eigenvalues ^from which 
the temperature multiplied by entropy is subtracted and then they are multiplied by the 
same projectors. If we consider the principal state of this quantum free energy, then we 
obtain all the results of the former thermodynamics. But we can consider not only the lower 
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level of this serial but also the series corresponding to the velocities (currents). Then we 
obtain answers to the question of how the metastable states vary with temperature. 

We shall show that there exist serials that contain superfluid and superconducting series 
at any temperature. However, these series necessarily do not contain the principal state, i.e., 
there exist serials whose lower level is lower. By idea, the nature must choose the latter. 
On the other hand, by the argumentation of the great Dirac, "it were astonishing if the 
nature did not use this possibility." I add "somewhere", in nuclei, in stars, etc. in one word, 
somewhere. 

Is it possible to create such serials "by hand"? We have already created by hand 
metastable series, higher than the principal state, by letting currents flow. Is it possible, by 
using resonance, also to "pump" high-temperature serials? 

1.4 Necessary conditions for metastable states 

If a large parameter is given as some concrete number, it is sometimes not clear whether this 
parameter is "large" or "not very large." Say, if the number of particles is 100, this number 
must be a sufficiently large parameter, but the logarithm of this large parameter is already 
not too large. In this case, first of all, we would like to know how to determine a metastable 
state. I introduce, at least, necessary conditions for a metastable state. 

To this end, I first consider a very simple example, the one-dimensional Schrodinger 
equation whose potential well has two hollows, i.e., there are two minima: one minimum is 
global, i.e., it corresponds to the principal state, and the other minimum is local. In the case 
of the classical asymptotics, it determines a metastable state. What does it mean? 

If the system is at this lowest level corresponding to the local minimum or to the smaller 
hollow, i.e., in other words, if the eigenfunction is concentrated near the local minimum and 
rapidly decreases everywhere else, then the following question arises: what happens if the 
system if perturbed? In other words, what is the probability of the transition ,Trom this 
state to the principal state, i.e., to the state that corresponds to the global minimum and 
whose eigenfunction is concentrated near this minimum and decreases sufficiently rapidly 
everywhere else. 

If we have a semiclassical approximation, then the tails of these eigenfunctions decrease 
very rapidly. Hence if we consider the matrix element obtained as a result of perturbation 
by some operator of multiplication by an external potential, then we see that this matrix 
element is very small, since the product of functions is very small. In other words, I take 
the following matrix element: ipo is the eigenfunction corresponding to the principal state 
and if) i is the first eigenfunction corresponding to the other hollow. V is the potential by 
which we perturb the system. If the matrix element (ip , V'tpl) is small, then this means that 
the transition probability is small. Hence this state is metastable, i.e., it is very difficult to 
"fall" to the lower level, while transitions to the higher levels whose supports are also near 
the minimum of this smaller hollow are quite possible. 

However, since we say that this metastable state is, to some extent, a model of super- 
fluidity and superconductivity, it is very important for us that the energy did not become 
less. If the energy increases, then the larger the better. (If the velocity increases due to 
this perturbation, this is also very good.) So, we perturbed the problem as if by friction, 
as Landau would say, and nevertheless, the energy did not decrease, since the transition to 
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lower levels for such an asymptotic problem is practically forbidden. 

However, if we take an operator not from the class I described but from a more com- 
plicated class, e.g., from the class where the translation operator is also involved, then the 
matrix element is not at all small. Hence the first problem is to specify the operators by 
which we perturb our problem. 

The operators by which we perturb the problem are operators that preserve the support 
of a function, i.e., in other words, if a function ip is sufficiently smooth and equal to zero 
outside some interval (thus is interval is called just the support of this function) and we 
apply the multiplication operator or the differentiation operator of any order, then these 
operators do not take the function outside its support. We usually perturb our problem by 
such operators. What does it mean? 

This means that, prior saying that we have a metastable state, we must define the class 
of operators with respect to which this state is metastable and for which the transition 
operators are small. Natural operators of this type are just the potential forces. Further, 
the question arises: perhaps, this is not a semiclassical approximation, the parameter h is 
not small, perhaps, the barrier is not very high, what happens then? Can a metastable state 
exist in this situation and how long? In this section we shall try to answer this question. 

The first question: is such a state possible at all or not? The second question: how small 
are the matrix elements of the transition? Thus the first question is reduced to the following 
one: does the minimum exist, or not, and of what is this minimum? 

We spoke about the minimum of a potential well. However, about what minimum must 
we speak in the quantum problem itself, without a small parameter? 

In the quantum problem without a small parameter we can reason in the following way: 
we consider a set of perturbing operators with support on the interval [a, b], in other words, 
we take the operators of multiplication by smooth finite functions that differ from zero on the 
interval [a, b}. Then we act by the set of such operators on the entire space L 2 . We obtain the 
set of values of the result of this action. We study whether the original Hamiltonian attains 
its local minimum on the functions from this set. Namely, this minimum is understood as the 
true local minimum with positive second-order derivatives, i.e., the minimum does not lie on 
the "edge" of the interval but lies somewhere in the middle. If the interval [a, b] includes the 
minimum point of our potential, i.e., the smaller hollow, then the local minimum necessarily 
exists for a sufficiently small H. Such a minimum can also exist if H is not very small, but just 
the existence of such a minimum answers our question of whether there exist a metastable 
state. 

If there is no minimum, we cannot speak about a metastable state. If there is no second 
hollow, then, apparently, the Hamiltonian does not have a minimum, at least for sufficiently 
small h. 

At all events, we, first of all, must answer the question of whether the quantum problem 
with classical Hamiltonian, which is not related to the classics, has a local minimum. If 
such a local minimum exists, then there exists an eigenvalue, which is the nearest to this 
minimum, and the corresponding eigenfunction determines some metastable state of the 
original operator. In this case a small perturbation of the length of the interval does not 
lead to the transition to another value. 

The value of the dimensionless variable Voa 2 m/h 2 at which the local (not principal) mini- 
mum disappears will be called the critical value. Here m is the mass, % is the Planck constant, 
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and the potential has the form Vov(x/ a), where Vq is a constant, a is the characteristic length, 
and v(y) is the potential well with two hollows mentioned above. 

The question of how long this metastable state lives is still open, hence I only say that 
this condition is necessary for the existence of a metastable state. 

1.5 Necessary conditions for the existence of superfluidity and 
superconductivity 

Obviously, a metastable state does not depend on the representation. However, we relate 
the problem of finding a minimum to some asymptotic problem, where the representation 
in which we seek the minimum (namely, the /c-representation) is chosen in some natural 
way. Similarly, we act in the many-body problem, in the situation in which the number of 
particles is not very large. We first determine the representation for which these metastable 
states are natural in our asymptotic problem. 

Hence in the many-body problem we choose the representation in which, as N — > oo 
we already have series with metastable lower levels. The most convenient way is to deal 
with the representation of occupation numbers, which is introduced in the supersecondary 
quantization, and hence, in the variational method, which will be considered later. 

These occupation numbers rij attain integer nonnegative values for j — 0, 1, 2, , . . .. 

Let us consider the set of functions of occupation numbers A(no,ni,ri2, ■ ■ .) such that 

\A(n ,n 1 ,n 2 , . . .)| 2 < oo 

M 

and A(n , rii, n 2 , . . .) = if N ^ n + n\ + n 2 + . . .. The subset of functions of occupation 
numbers such that A(n ,ni,n 2 , . . .) = for aN < n + ni + n 2 + ■ ■ ■ will be the characteristic 
class of functions corresponding to the number < a < 1. 

To determine a metastable state in the problem of N bodies, we seek the minimum of 
the Hamiltonian in the representation of occupation numbers on the set of functions from 
the characteristic class. By analogy with the preceding section, we can find the critical value 
N CT for which the metastable state disappears. Superfluidity and superconductivity in the 
system are possible for iV larger than the critical values. Thus, iV cr is the lower bound for 
N for which superfluidity and superconductivity are possible at least for a short time. 

1.6 Variational methods and variational principles 

There exist many variational methods, for example, the variational Riesz method, the varia- 
tional Galerkin method, the variational Tamm-Dankov method. The method of finding the 
solution is, so to say, successive, i.e., there a priori exists a sequence that converges to the 
solution. In the Tamm-Dankov method this convergence is not proved, but each successive 
approximation must be, in principle, better than the preceding approximation. 

The situation with variational methods is somewhat different. The variational methods 
are formulated so that they directly lead to some equation and are postulates that follow 
from the equations. There also exist variational principles that are derived from the original 
variational principle which is considered on a narrower function class. They lead to some 
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new equations. The relation between these new equations and the original "true" equations 
is determined experimentally. 

Suppose that some variational method, say, the Hartree-Fock method is used. The 
results obtained are compared with experimental results, and the coincidence is good. In 
the variational principles such as the Hartree-Fock principle, the variational Bogolyubov 
principle, the variational Schriffer principle, we try to guess the form of the solution by 
substituting this solution and obtaining some new equation. So, these principles lead to some 
new equations and, instead of the original equations, we solve some other equations obtained 
from the variational principle. We could determine the relation between the variational 
principle and the equations if we had a sequence (as in the Riesz method, in the Galerkin 
method, and in some other methods) of test functions with the help of which we could step- 
by-step approach the solution, but not a single test function with the help of which we need 
to verify the variational principle. 

Just in this case, the supersecondary quantization, which I propose, naturally leads to 
a variational method, i.e., to the method of successive approximations. To some extent 
this ideology is close to that of the Tamm-Dankov method. The variational Bogolyubov 
and Hartree-Fock principles constitute only the first approximation. For some problems 
this yields exact asymptotic methods, as was shown by Bogolyubov for the Cooper-Bardeen 
model. He showed that this is an exact asymptotics as iV — > oo. It is of interest that the 
variational method proposed here leads to better results if the number of particles is not 
very large, just where the asymptotic theory cannot be applied. 

If the number of particles is 104, as in traps, then this case cannot be an asymptotics 
as iV — > oo, since the logarithm, which tends to infinity as N — > oo, is here equal to 4. 
If we have 10 23 , just as in the weakly nonideal Bose gas, then the logarithm is already a 
sufficiently large number, and this case can be an asymptotics for large N. In this case it is 
not so convenient to use the variational method. 



2 Supersecondary quantization in the boson case 

Following [3], we introduce the averaging of operators in the more general case of a large 
number of two types of clusters. Let k > 1 and N\, . . . , N*. > be integers. Consider the 
Hilbert space of functions of the form 

,w 1-1 1 -1 2 2 -2 2 2-2 

^{ x li Jl'i • • • i x Nii J Nil x li x 2iJli ■ ■ ■ i X 2N 2 -1' X 2N 2 '3n 2 i ■ ■ ■ 

x li ■ ■ ■ i x ki3\ }■■■'■> x kN k -k+lt ■ ■ ■ i x kN k iJN k )i (1) 

where x l p[ G T 3 , T 3 is the three-dimensional torus with side length L, j l m = 1, . . . , oo for all 
I = 1, . . . , k and p\ = 1, . . . , INi, qi = 1, . . . , Ni. moreover, the functions (§) are symmetric 
with respect to the permutations of any pairs (x l pl _ l+1 , . . . , x l pl , j l p ) and (x l ql _ l+1 , . . . , x l ql , j l q ) of 
sets of variables, where I = 1, . . . , k and p, q = 1, . . . , N t . 
The inner product has the form 



n n 

= E ...£••• E ... dx\...dx\ l ...dx\...dx 



kN k X 
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(x\, j\; . . . ; x x Nl , ; . . . ; xj, . . . , xjjl, j\ ; . . . ; x kNk _ k+1: • • • ■> x kN k i3N k ) 

A ^l x H Jl i • • • ) X N 1 ^JN 1 J ■ ■ ■ i X D • • • ) x ki J I i • • • J x kN k -k+li • • • i • L kN k iJN k ) 

We denote this space by J r N 1 ,...,N k - Let us define the boson Fock space F§ N . 
Definition. The space J-^ N is the subspace of J r N 1 ,...,N k formed by the elements 



symmetric with respect to arbitrary permutations x l v <-> x™ for l,m = 1, . . . , k and p 

dB 

N u ...,N k - 



1, . . . , INi, q = 1, . . . , mN m . The projection on this subspace will be denoted by Pjf 



D N 1: ...,N k '■ F Nl ,...,N k -> ^ lv ..,JV fc - 



Now consider the Hilbert space whose elements are infinite sequences 

* = {*iv 1) ...,ivj, e ^JV!,...,^, JVi, . . . , JV fc = 0, 1, . . . , 

oo oo 

u*ir= e ••• e ii^jvi jv*ir < oo- 

JVi=0 JV fc =0 

This space is the infinite direct sum 

? = F Nl ,...,N k 
N u ...,N k 

and is a special case of the Fock spaces [4]. 

In this space we introduce the creation operators b'l(x 1 , . . . ,xi,j) and the annihilation 
operators bf(x 1 , . . . , xi,j), I — 1, . . . , k. 

Let PN u ...,N k be the projection on the corresponding component of the direct sum, 

P N u ...,N k ■ F -> FN u ...,N k , PN lt ...,N k y = *JVi JV fc , 

and let iN lt ...,N k be the operator of embedding of J-Ni,...,N k in 

*JVi,...,JV fc : F Nl ,...,N k — >■ ^ , iN!,...,N k <P = {^JViMi • • • • • 5 Nk M k <p}, 

where is the Kronecker delta. 

We introduce the notion of an averaged operator. Let {^4tv} be a sequence of operators 
in the subspaces of L 2 ( T 3 x . . . x T 3 ) formed by elements ip(xi, . . . ,xn) symmetric with 

n times 

respect to the permutations of the variables Xj and x k , Xj € T 3 . To these operators, we 

assign the operators A Nlj _ jNk , Ni + 2N 2 + . . . + kN k = N defined on the subspaces F§ N 
and acting as the identity operators with respect to the discrete variables j l . 

Definition. The averaged operator A is the operator 

oo oo ~ 

A = E • • • E C^l + • • • + kN kV- i N 1 ,...,N k A Nl+ ... +kNk P Nu _^ Nk P Nu ... tNk 
N 1= N k =0 
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in the Fock space T . 

Remark. The operator i Nl ^„ >Nk AN 1+ ... +k N k PN- i ,...,N k PN 1 ,...,N k represented as 

1 OO OO OO OO 

*rj /v~7 E ■•• E ••• E •• E / •••/ ^...dx^.-.d^...^ 

J2 _2 J2 

x ^1 (^1) Jl) • • • &l \ x Nii3Ni) ■ ■ ■ \ ( X li ■ ■ ■ ' X ki 3 k) ■ ■ ■ 

2 

x b k (x kNk _ k+1 , . . . ,x kNk ,j Nk )A Nl+ ^ +kNk 

l i 
x Symm { 6" (ar} , j{) . . . b± {x]f , j x Nl ) . . . 

11 k k 

X 1--- X N 1 ■■■ X l--- X kN k 

1 1 

x fyc • • • i x \i3k) ■ ■ ■ b k ( x kN k -k+i> ■ ■ ■ ? x kN k 5 3% k ) } 
x exp ( - ^ / dz bf (z,j) \ (z,j) - ... 

oo „ „ 2 1 

_ E / ••• / dz i--- dz k bt (z!,...,z k ,j)b k (z!,...,z k ,jf), (2) 

where Symm a . l _ is the operator of symmetrization with respect to the variables Xi, . . . ,xn 
and numbers over operators indicate the order of action of these operators [5]. The expression 
(H) was obtained in [6] for the special case k = 2 and can be proved in a similar way for 
arbitrary k. 

According to this remark, the following functional will be called the symbol of the averaged 
operator A: 



A{bl{x{, j 1 ), hixlJ 1 ), b* k (x%, x k k ,j k ), b k (x1, x k k ,f)) 



E---E 



(Ni + 2N 2 + . . . + kN, 



Ni=0 N k =0 ^!.iV 2 iv fe . 



X 



OO OO OO p » 

x E • • • E • • • E • • • E j ■■■ j dx\... dx l Nl ...dx k ... dx k kNk x 

x ^i( a; i? Ji) • • • ^i ( x jvi ' i_Vi ) • • • b k (x k , • • • , x k ,j k ) ■ ■ ■ b k {x kNk _ k+1 , . . . , x kNk ,j Nk ) x 
x ^v 1+ ...+fcv fc Symm {feifcLji) . . , ) . . . x 

1 1 k k 



x l— x N 1 — x l— x kiM h 



xb k (x k , . . . , x k ,j k ) . . . b k (x kNk _ k+1 , . . . , x kNk ,jx k )} 

OO „ 

x ex P ( - E / dzb *i(z,j)bi(z,j) - 

OO „ - 

~E / ••• / dz 1 ...dz k b* k (z 1 ,...,z k J)b k (z l ,...,z k J)y (3) 



9 



where &/(•, . . . , € L 2 (T 3 ), I = 1, . . . , k, for all = 0, 1, . . . and the inequality 
E / •• / dx 1 ...dx l b* l {x 1 ,...,xi,j)b l (x 1 ,...,x h j) < oo 

is valid. To the sequence of operators {An}, there corresponds a secondary quantized oper- 
ator A in the Fock space Ji B - 
The following assertion holds. 

Lemma. The symbol of the averaged operator A can be expressed as follows: 

^ oo - 

= Sp {Ap) exp ( - ^ / dztiftzjfaizj) - 

-••• _ E I ■■ I dz 1 ...dz k b* k (z 1 ,...,z k ,j)b k (z 1 ,...,z k ,j)), (4) 

3=1 

where p is an operator in He and has the form 

oo oo i k oo „ „ 

p= E ••• E jrr~jri IT(E /••• dx[d y [...dx\dy\x 

x6i*(2/[, • • • $~ (yi) • • • (y!)) ! ex P ( -y ^ ^ + (*) ^" (*)) • 

This lemma follows from a lemma proved in [3]. 

As was shown in the Introduction, it is impossible to find how metastable states depend 
on the temperature, if the quantization of free energy is not introduced. 

In the simples case of microcanonical distribution, we have the following expression for 
free energy F: 

F = X-6\nn x (5) 

where 9 is the temperature, A is the energy level, ri\ is its multiplicity, and S = lnn x is 
entropy by definition. Let dE\ be the spectral family of a self-adjoint energy operator H 
with discrete spectrum. Then, as is known, we have 



H = J XdE x . 

It is natural to define the quantum free energy as 

F = J (A - 9\nn x )dE x . (6) 

As a rule, in statistical physics we deal with the canonical distribution. The passage from (5) 
to the canonical distribution is not trivial and here can be carried out only in the simplest 
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cases. Similarly, the method for obtaining the values of quantum entropy, which we propose 
later, can be derived from (6) not in the general case in which we shall use it. This separate 
problem lies beyond the framework of the present work. 

Next, let Ti and E be the Hamiltonian and the identity operator, respectively, averaged 
by this method. By the lemma, the symbols of these operators, up to a factor, are equal to 
Spp and Sp(Hp), respectively. If p is treated as a statistical operator (a density matrix), 
then Spp is the norm of this operator and Sp(Hp) is the mean value of the Hamiltonian 
multiplied by Sp p. Furthermore, entropy corresponding to p is given by 

s = Sp ( ?ln (sf?)) (Sp?r1 ' 

By analogy with statistical physics, we refer to the functional 

S(bl(x{, j 1 ), hixlJ 1 ), . . . , b* k (x k , . . . , x k k ,j k ), b k (x k , . . . , x k ,j k )) = 
= Sp (pin ) ex P ( - £ / dzbi(z,j)b!(z,j) - ... - 

CO „ „ 

~S / ••• / dz i ■ ■ ■ dz kK( z i, ■ ■ ■ , z k, j)bk(zi, ■ ■ ■ , Zk, j)) (7) 
j= i J J 

as the symbol of averaged entropy. To this symbol, we assign the following operator in T: 

§ = RegS(bt (xlj 1 ),^ (x\,j 1 ),...,bt(x k ,...,x k ,j k ),i k (x k ,...,x k ,j k )), (8) 
where Reg stands for the regularization if necessary. 

Definition. Numbers A such that 

(7? + 6S)§ = AE$, $G $ ^ 0, (9) 

are called quantum eigenvalues of free energy for a system of bosons with Hamiltonian H at 
temperature 9 for k cluster types. 



3 Supersecondary quantization for fermions 

Now we introduce the averaging of operators and the quantum free energy for fermions. 
We assume that the fermions may have a spin variable s ranging in a discrete finite set E. 
Consider the Hilbert space Fn-l,...,^ of functions (1), where x l pi G T 3 x E. Throughout this 
section, we use the notation 

/ dx = £ / de, (io) 

where x = (£, s), £ G T 3 , s G E. Let us define the fermion space 

Definition. The space J 7 N 1 ,...,N k * s the subspace of J r N 1 ,...,N k formed by elements J r N 1 ,...,N k 
that are antisymmetric with respect to the permutations of arbitrary x l p and x™, l,m — 
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1, . . . , k, p — 1, . . . , INi, q = 1, . . . , mN m . The projection on this subspace will be denoted 

pF . r fpB 

r N u ...,N k ■ S Nl ,...,N k -> ^JVi,...,JV fc - 



Let us introduce the notion of an averaged operator in the fermion case. Let {A N } be 
a sequence of operators in the subspaces L 2 (T 3 xEx...xT 3 xS) formed by elements 

S v ' 

n times 

ip(xi, . . . ,xn) antisymmetric with respect to the permutations of the variables Xj and Xk, 

Xj G T 3 x S. To these operators, we assign the operators Ajv!,...,^, Ni+2N 2 + ■ . . + kN k = N, 
defined on the subspaces J~Ni,...,N k an d acting as the identity operators with respect to the 
discrete variables j 1 . 

Definition. The averaged operator A is the operator 

OO CO „_ 

A= E ••• E(^i + --- + ^)^ 1 ,..,Ar fc A iVl+ ... +jfci v fc P^ ) ... jJVfc P A r 1 ,..., A r fc (11) 

JV 1= N k =0 

in the Fock space T . 

Remark. The operator iNi,...,N k A.N 1 +...+kN k PN 1 ,...,N k ^ > N 1 ,...,N k can be written as 

OO CO CO CO „ . 

n W~\ £ ••■ S ••• E ••■ E /•••/ dx\...dx l Nl ...dx k 1 ...dx 

'l--" iV *- ji=l .,1=1 ,-*=! ,t r =1 



k 

^ ... ^ ^ J ■■■] 

J2 _2 ^2 

x bf (x 1 . . . (x Nl , j Nl ) . . . b k (x l , . . . , x k , j k ) 



X ^fc (^fcAffc-fc+l? • • • > x kN k i jN k )AN!+...+kN k 

x Asymm { b± (x\ bi (x 1 ^ , j 1 ^) . . . 

■ L 1--- L N 1 ■■" L l--- L kN k 

J. 1 

x b k (x k , . . . ,x k ,j k ) . . . b k (xl Nk _ k+v . . . , x kNk , J^)} 

CO r J J 

x exp ( - E / dz bf (z,j) \ (z,j) - ... 

-E / ••• dz 1 ...dz k bt (z 1 ,. . . ,z k ,j)b k (z u . . . ,z k ,jf), (12) 

where Asymnx^ is the operator of antisymmetrization with respect to the variables 
xi, . . . ,xn- This expression is the fermion counterpart of (2). 

The symbol of the averaged operator A in the fermion case is defined as the functional 

A{b\{x\, j 1 ), biixlj 1 ), b* k {x k , x k k ,j k ), b k (x k , ...,x k k , /)) = 
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~ ~ (N t + 2N 2 + . . . + kN k )\ 

2^ ■■■ 1^ /vj/Vo' 7VJ x 

AT 1= Af fc =0 JVi.JV 2 iVfc. 



Vfc = 

oo 00 



p 

x E • • • E • • • E • • • E J ■■■ J dx\... dx 1 ^ ...dx k ... dx k kNk x 

x ^ ) l( ;r l) Jl) • • • W (icjvi ) Jn-l ) • • • bk( x li ■ ■ ■ 5 x ki3k) ■ ■ ■ bk( x kN k -k+li • • • 1 x kN k i3N k ) x 

x^7V 1+ ...+fcJV fc Asymm {feiO^Ji) . . .&i(a^ ) . . . x 

• i i--- x iv 1 •■' x fcJV fc 

xb k (x k , ■ ■ ■ 1 X kiJk) ■ ■ ■ bk( x kN k -k+l-> ■ ■ ■ 1 X kN k 5 3N k ) } 

00 . 

x ex p ( - E / dzb l( z ,j) b i(z,j) - 
3=1 

~E / ••• / dz 1 ...dz k b* k (z 1 ,...,z k ,j)b k {z 1 ,...,z k ,j)y (13) 
where . . . , e L 2 (T 31 ), I — 1, . . . , k, j l — 0, 1, . . ., and the inequality 

OO i, „ 

E / • • / da;i-..da;i6r(a:i,...,xj,j)6j(xi,...,Xi,j) < 00 

is valid. To the sequence {^4jv} of operators, there corresponds a secondary quantized opera- 
tor A in the fermion Fock space Tip, which consists of sequences of antisymmetric functions 
and is an analog of the space TLb (see [4] for details). 
The following assertion holds for the fermion case. 

Lemma. The symbol of the averaged operator A is given by the formula 

A(bl(x{, j 1 ), hix],] 1 ), b* k (x k , . . . , x k k ,j k ), b k (x k , ...,x k k , /)) = 
= Sp (Apj exp ( - J dzbl(z,j)b!(z,j) - 

-...-E / ••• / dz 1 ...dz k b* k (z 1 ,...,z k ,j)b k (z 1 ,...,z k ,j)), (14) 
j= i J J 

where p is the operator in Tip given by 

00 00 00 00 00 00 

p= E ••• E N , 7Y 1 E • • • E ■ E >< 

jv 1= o Jv fc =o Jv i---- iV fc- i i =0 ^=0 j*=0 j^=0 



x J ...J dx\dy\ . . . dx^dy 1 ^ . . . dx\dy\ . . . dx k kNk dy k kN J) + {x\) . . . ^ + (^j Vl ) x 
x...^(x k )...^(x k kNk )P x 

x bi(x\,jl) . . . bi(x]^ i: jlfj . . . b k (x k , . . . , x k ,j k ) . . . b k (x kNk _ k+1 , . . . , x k N k i3N k ) x 
xbUylJl) . . . b\(y l Nv3 l Ni ) ...b* k (y k ,..., y k ,j k ) . . . b* k {y k kNk _ k+1 , y k kNk ,j k Nk ) x 
x$-(y k kNk ) . ..r{y k ) . ..^(Vk) ■ ■■#-(*!); (15) 
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here P is the projection on the vacuum vector in Tip- 



In the fermion case, we define the symbol of averaged entropy as the functional 
SibKxlj 1 ), hixlj 1 ), b* k (x k , x k k ,j k ), b k (x k , x k k ,j k )) = 
= Sp |pln (g^^ ex P ( - £/ dzb *i(zJ)bi(z,j) 

-J2 / ••• / dz 1 ...dz k bl(zi,...,Zk,j)bk(z 1 ,...,z k ,jj) (16) 
j= i J J 

To this symbol, we assign the operator 

S = Re g S(if (xlj 1 ),^ (xlj 1 ),...,^ (x k ,...,x k ,j k ),b\ (x k , . . . ,x k , j k )) (17) 
in the space T . 

Definition The numbers A such that 

(H + 6S)§ = \E<Z>, $ G T, $ ^ 0, (18) 

are called quantum eigenvalues of free energy for a system of fermions with Hamiltonian H 
at temperature 9 for k cluster types. 

4 Asymptotic series as N — » oo for a fermion system 

Further, we assume that the secondary quantized fermion Hamiltonian in the space Tip has 
the form 

H = J! dxdy^ + (x)T(x,y)^-(y) 

+IIIII dxdx'dydy 1 ' $ + {x)$ + {x')V{x,x' ,y,y')$-{y')$-{y), (19) 

where, just as everywhere, x,x',y,y' = ((,s), ( G T 3 , s G S, and the functions T(x,y), 
V(x, x', y, y') have the properties 

T(x, y) = T*(y, x), V(x, x' , y, y') = V*(y, y', x, x'). 

Let us consider the super secondary quantized equations for free energy of fermions in another 
representation. Suppose that k is the number of clusters of the first type, M is the number 
of clusters of the second type, and there are no clusters of any other type. We choose some 
functions V(x,y) G L 2 ((T 3 x £) x (T 3 x £)), (j)(x,j) G L 2 (T 3 x £) x l 2 such that 

q?(x,y) = -V(y,x), 

fc = Sp^G* - 1 + (l + (1 - G*) _1 ^(l - G^R+y 1 ^ , 

l - Sp^(l + (1 - G^Ril - G^R+y 1 ^ - G')- 1 ^! - G)- 1 ^, (20) 



M 

2 
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where G, R are operators in L 2 (T 3 x E) of the form 



G = (1 + A)~ yA f + B + (1 + A)~ L B(1 + A f ) [1 + B + (1 + A)~'B{1 + A 
R = -(1 + A^Bil + A*) _1 (l + £+(1 + ^-^(l + A*)- 1 )" 1 , 
and A, £> are the following operators in L 2 (T 3 x E): 



>u-i 



(21) 



(-Au)(x) = / d yp( x J) ( P*{yJ) u (y)^ ( Bu )( x ) = / dy^(x,y)u(y), 
j=i J J 

where the superscript + stands for Hermitian conjugation and t for transposition. We choose 
an arbitrary complete orthonormal system of functions ty a (x,y), a = 1, 2, . . ., in a subspace 
of the space L 2 ((T 3 x E) x (T 3 x E)) whose elements are functions orthogonal to ^(x, y), and 
an arbitrary complete orthonormal system of functions ip c (x,j), c = 1, 2, ... in the subspace 
of the space L 2 (T 3 x E) x l 2 consisting of functions orthogonal to ip(x,j). The following 
vector system corresponding to these systems is complete and orthonormal in the subspace 
J~kM of the space T 



{n},{m} 



= n 



o=l V^a! 



dx 1 dx 2 ^ a (xi,x 2 )b 2 (xi,x 2 ) 



oo 



1 



c= i v«v y j=1 

1 



X 



X 



I d y<Pc(y,j)bt(y,j) ] x 



dxdx'*5/(x, x')b 2 (x, x) 



I oo \ ^ 

- 1 E / dyW,f)bt(y',f) *o, 

d n iJn ± \ \r=i J J 



(22) 



where $ is the vacuum vector of the space T\ {n}, {m} are sets of integers n a > 0, 
a = 1, 2, . . ., and m c > 0, c = 1, 2, . . ., such that 



(2=1 C=l 

and the following notation is used: 

iV 2 = M-f> a , iV 1 = A;-f> c , 

and 



a=l 



c=l 



C = \j II dxdX ' l^^'^')' 2 ' D = \^J2f dx \f( x J 
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We introduce the Fock space M. generated by the vacuum vector Y and the following boson 
creation and annihilation operators of two kinds: (3^, a = 1, 2, . . ., df, c = 1, 2, . . .. To each 
vector (22) we assign the following vector of the space M: 

oo -i OO -I 

%w = n -iTnCtiY* n 7^r^) mcF o- ( 23 ) 

o=l V u a- c=i V ''*<:■ 

Taking into account the relation between (22) and (23), we can write the equation for 
quantum values of free energy of fermions in the space M. as 

F e Y = XEY, Y e M, Y^O, (24) 

where Fg, E are the following operators in M.: 

oo 2 1 oo 2 1 



F e = exp -EA + A~-E rf X E E ^},{™}-( W + ^) $ M,M 

V i=l P=l / {n},{m} {n'},{m'} V 

° (fr)"-(fr) Bi ) ^ (4) mc (t) m ' c ) 



a X 



X 



oo J J oo 2 1 

^ = exp I - E A + A~ - E d+d~ I E E ( $ W,{m}, 



X 



i=l p=l / {n},{m} {n'},{m'} 



x 



i t , t: - n 



a =i \Jn a \n' a \ c =i yjm c \m' c \ 
We assume that equation (24) has the solutions 

Y = E l{n},{m}Y{n},{m}, (25) 
{n},{m} 

such that in the limit as N — > oo, L — > oo, iV/L 3 — > const we have 

E{n},{m}(E^l^a)|7W,{m}| 2 _ . s E{n},{m} (Egll m c ) |7{n},{m} | 2 v 

El'v I 2 ^V- 1 -/' I |2 *-A-U" 

{n},{m} 1 7{n},{m}| 2^{n},{m} I l{n},{m}\ 

To such solutions of equation (24) there correspond the quantum value of free energy 
K^( x ^y)^( x J)) = II dxdyT(x,y)G(x,y) + 
+ 2 



lIIII dxdx ' d v d y' v ^ x ' ^y^'){ G ^yW ^y') 

G(x, y')G(x', y) + R*(x', x)R(y', y)) + 



+^S P /( ° ^ )+0(l), (26) 
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where /(£), (6R has the form 



/(0 = 5(5 + + + 5(5 - 1*4 " 0, 

and the functions G(x,y), R(x,y) are the kernels of the operators (21): 

(Gu)(x) = / dyG(x,y)u(y), (Ru)(x) = / dyR(x,y)u(y). 



The vector $ e corresponding to (p5|) we denote by ^(^(x, y), <p(x,j)). For two distinct 
pairs of functions *$?(x, y), ip(x,j) and *ff'(x, y), ip'(x,j) satisfying conditions (p0|) , the follow- 
ing matrix element of an arbitrary supersecondary quantized operator A, which is bounded 
in N: _ 

$(*'(x, y), (f'(x,j)), A${V(x, y), <p(x,j)) 



W*'(x,y),if/(x,j)),$(*(x,y),(p(x,j))) 

is exponentially small in N. 

Therefore, in the case in which the values X(^(x,y),<f(x,j)) and A(\E ,/ (x, y), <p'(x, j)) 
corresponding to these function pairs do not coincide, the vectors $(^(a;, y), (p(x, j)) and 
$(\P'(x, y), ip'(x, j)) belong to different series. If the resonance occurs, i.e., the value (|26|) is 
the same for *&(x, y), ip(x,j) and ty'(x,y), tp'(x,j), then the series is formed by linear combi- 
nations of the vectors $(\l/(x, y), <p(x,j)) and $(\l/'(x, y), ip'(x,j)). From all series of solutions 
of the equation for quantum values of free energy corresponding to some chosen k and M, 



we choose the series with the least value (|26|) . The functions ^(x, y), tp(x,j) corresponding 



to this series are determined from the equation of self-consistent type: 

G-\ R + \ 1 , , ( 1 ( T V 



R -G % + i J 2 \ 29 



where the operators T, V are expressed in terms of G, R as follows: 
T = X - fi + (2// - uj){1 - G l )- l R(l - G)- 1 J R+(1 - G 1 )' 1 



x 



(i + i?(i-G)- 1 i?+(i-G*)- 1 



-2 



V = Z + (2fi-uj)(l-G t )- 1 R(l-G)- 1 (l + R + (l-G t )- 1 R(l-G)- 1 ) , (28) 

and the operators X, Z in the space L 2 (T 3 x l 2 ) have the form 

{Xu){x) = J dyT(x,y)u{y) + JJJ dx'dydy'(V(x,x',y,y') -V{x,x',y',y))G(x',y')u(y), 
(Zu)(x) = III dx'dydy'V(x',x,y,y')R(y',y)u(x'). (29) 



The numbers /i, uo in ( |2"8"D are determined by conditions (^). The minimal quantum value of 
free energy ([26|) for some chosen k and M we denote by \u,m- Let us consider the minimum 
of Xk,M with respect to k and M provided k + 2M = N = const). We denote this minimum 
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by Ajv. The functions *&(x, y), tp(x,j) corresponding to this minimum can be found from the 
equation 



G-\ R+ \ = _ l _ ih L Xt -^ f 
R + | ) 2 \26 \ Z -X + n 



(30) 

where the operators X, Z are the same as (|2~9| ) and the number [i is determined by the 
condition k + 2M = iV, where k and M have the form (|2~0"D . Equation (|30|) is well known. It 
coincides with the Bardeen-Cooper-Schriffer-Bogolyubov equation in the theory of super- 
conductivity. The solutions of this equation with different 9 determine the functions k(6), 
M{9) by formulas fl2"0|) and the function \n{0) by formula (f26|). The temperature 9 C at which 



M{9 C ) = is called critical. At this temperature the heat capacity 

C = -6\" N {6) (31) 

has a discontinuity. Since resonances can occur for Xk,M = ^k',M' with k ^ k! and M ^ M', 
we can replace the functions k(6), M{9) by some other ki(8), Mi(8) such that Mi(#i) = for 
6i > 6 C and thus to increase the critical temperature. To calculate the heat capacity corre- 
sponding to a temperature series distinct from the Bardeen-Cooper-Schriffer-Bogolyubov- 
series, we need to replace \n{9) by Xk(e),M(e) i n formula (|31|). 

5 Asymptotic series as N —> oo for a boson series 

In a similar way, we can find series of quantum values of free energy for bosons in the case 
in which the number of cluster types is 2. 

Here we present an asymptotics for the case in which the condensate can exist. 

The principal state of the series is expressed in the same way as for fermions via the 
vectors of some system of the form (21). For the sets {n}, {in} the same condition as that 
for fermions holds in the thermodynamical limit. Let us study the series such that 



1^0(p) e ^), 0(p) = 0(-p), ^j) = J^M e ^ 



p 



where p are three-dimensional vectors of the form 27r(/i, l 2 , h)/L, li, / 2 , / 3 G Z, J2 P stands for 
the summation over all such vectors, and pj is a single- valued mapping of the set j = 1,2,... 
on the set of such vectors. Moreover, 0(g), a(q) are continuous functions of the variable 
q G R 3 such that 0(g) = 0(— g), a(g) = a*(g), < a(g) < 1. If the condensate consists of 
particles with zero momentum, then the functions 0(g), a(g) satisfy the inequality 

(1 - a(g))(l - a(-g)) - 0*(g)0(g) > 0. (33) 

Note that we have the equality only for g = 0. In this case the asymptotics of the ratio of 
the quantum value of free energy to the volume in the thermodynamical limit has the form 

A 1 f, h 2 q 2 U a 3 n v(0) U a 3 n f 

* 2^ G{q) + 2 + W J dq {V{0) + < a ^ G ^ + 
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/ dqd{aqmq) + R * iq)) + J I dqdl im + ^ q ~ 0))g(?)g(0 + 

+ / / WtWq ~ l))R*(q)R(l) + 
6 r 

+ (2^p J dq(n(q)\nn(q)-(l + n(q))\n(l + n(q))), (34) 
where n is the condensate density, 

N is the number of particles, and just as for fermions we have N = k + 2M. The functions 
G(q), R(q), n(q) in formula (34) are expressed via a(q), 0(g) as follows 

a(g)(l-a(-g))+0*(g)0(g) 



G(q) 
R(q) 



[l-a(q))(l-a(-q))-<t>*(q)<f>(qy 

m 

;i-a(g))(l-a(-g))-0*(g)0(g)' 



1 



n(q) = g y/0- + G(q) + G{-q)f - 4R*(q)R(q) - -(1 - G(q) + G(-q)). (35) 

Expression (34) corresponds to the quantum free energy of the principal state of some series 
if the functions 0(g), a(q) satisfy the relations 

-a fn\ , L3 fj a{q) + a{-q) - 2a(q)a(-q) 



k = L3no ° (0) + 2(2^ I dq 



2{2ixf ./ * (1 - a (g))(l - a(-g)) - 0*(g)0(g) ' 

M = ^(1 - a(0)) + / rfg 7 , ^ — — (36) 

2 1 W; 2(2vr)3 7 y (i- a ( g ))(i_ a (-g))-0*(g)0(g) ^ ' 

and (34) attains the minimum on these functions provided that conditions (36) are satisfied. 
One can easily see that if n ^ then 

5 A 



U a s n v(aq) _^ Q 



, =0 2(27r)3(l-a(g))(l-a(-g)) 

This implies that for any function a(p) there always exists a function 0(p) such that the 
value of free energy (34) determined by the functions a(p), 0(p) is less than the value of free 
energy for a(p), <p{p) = 0. Although the nature can "set" the system on any metastable 
state, the physicists always assume that the state with the least energy is realized. Hence 
here we do not consider boson systems with clusters only of the first type and we always 
assume that there exists at least two types of clusters. 

6 Variational method for the supersecondary quan- 
tized problem 

We consider a variational method which allows us to construct successive approximates to 
the solutions of the supersecondary quantized equations for the quantum free energy both 
in the boson and fermion cases. We study the fermion case in detail. 
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Suppose that the number of clusters of the first type is k, the number of clusters of the 
second type is M, and there are no clusters of any other type. We study some complete 
system of functions ty a (x,y), a = 0, 1, . . ., in the space L 2 ((T 3 x E) x (T 3 x E)) and some 
complete system of functions ip c (x,j), c = 0, 1, . . ., in the space L 2 (T 3 x E) x l 2 . In this 
section we choose function systems no in the same way as in Sect. 3. Below we explain how 
we choose these function systems. 

To these functions, by formula (22) we assign a vector system that is complete in the 
subspace Tk,M of the space T . We assume that the numbers k and M are sufficiently large 

so that it is possible to regulate the entropy operator S. Let us consider the sequence 



$ b , [n + es $ 6 

\ b = ± )„ \ 6 = 0,1,..., (37) 



where G T is a vector of the form (22) corresponding to the set of numbers n a = M5 a b 
m c = k5 c b, and S a b is the Kronecker delta. We choose function systems ^ a (x, y), a = 0, 1, . . . 
and (p c (x,j), c = 0, 1, . . ., so that the following conditions hold: 

6Xb o, -^- = o ( 38 ; 



5V b (x,y) Sip b (x,j) 



and, in addition, 



min \ b = A . 

b 



Conditions ( |38D are Hartree-Fock-Bogolyubov-Bardeen-Cooper-Schriffer type equations for 
the functions ^(x,?/), (pb(x,j) = 0. 

For example, for ty(x, y) = at zero temperature, we can write the equation for the 
function ip b (x,j) as a system of Hartree-Fock type equations for the functions Uj(x) = 
<fb(x,j). 

In the equation for the quantum free energy we expand the vector $ with respect to the 
vector system (22): 

$ = ^{n}{m}^{n},{m}, (39) 

{n},{m} 

where 

E 

{n},{m} 

stands for the summation over all possible sets of {n}, {m}. Taking into account (|39]), we 
can rewrite the equation for free energy as the linear system 

E F{n'}{m'},{n}{m}A{n}{rn} = ^ E ^{n'}{m'},{n}{m}^{n}{ro}> (40) 
{n},{m} {n},{m} 

where the coefficients -F{n'}{m'},{n}{m}, -E{n'}{m'},{n}{m} are expressed via the supersecondary 

quantized operators Ti., S, E and the vectors (22). Equation (flCf ) is always the exact rep- 
resentation for the equation of quantum free energy. Approximate solutions of the equation 
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for free energy can be obtained from the linear system 




(41) 



where 0</c</c,0<M<M are integers and 



k,M 



E 



{n},{m} 



stands for summation over the sets of numbers {n}, {m} such that 



oo oo 



(42) 



a=l c=l 



For k = 0, M = 0, relation (41) implies the upper bound for free energy, which is equal to A 
(|37D . For k > 0, M > 0, this estimate can be improved and, for k — > k, M — > M, coincides 
with the exact value of the quantum free energy. 

7 Metastable states and the case of many clusters 

The same procedure can be used for any number of clusters both in the fermion and boson 
cases. Namely, we construct the sypersecondary quantization for this set of clusters and 
set the creation and annihilation operators equal to functions (i.e., we assume that they 
commute, and hence we take the symbols of these operators; the physicists say: " we assume 
that the creation and annihilation operators are '-numbers). 

Next, according to this nonlinear problem, we construct a set of independent functions 
complete in the corresponding space, which, in turn, determines the passage to the occupa- 
tion numbers. Then we employ the ideology of the Tamm-Dankov method. 

The set of such functions is constructed just as above, namely, the vector system (22) is 
naturally generalized to the case of an arbitrary number of clusters and the ratios of matrix 
elements similar to (37) are studied for this system. 

Expression (37) is the upper bound for the least quantum value of free energy. Hence 
it is natural to choose the function system so that this expression for each function attain 
its extremum value. Previously, this was performed in such a way that equation (38) is an 
extremum equation. For an arbitrary number of clusters, the approximate solutions of the 
equation for free energy can be constructed by using a variational method similar to (41). 

Hence, in some sense, this method is similar to the Tamm-Dankov method, since inequal- 
ity (42) is used. This variational method is especially effective for systems of third quantized 
fields and automatically can be used for such systems. 

The supersecondary quantization method can be used in quantum electrodynamics and 
quantum chromodynamics, when the "charge" is preserved rather than the number of par- 
ticles. Such an example is studied earlier. However, in this case we encounter the renormal- 
ization problem and must construct our variational method following the method proposed 
by Dyson [8], i.e., following the "new Tamm-Dankov" method. 
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To find metastable states corresponding to superfluidity and superconductivity, we per- 
form constructions similar to those considered in the Introduction. Let < a < 1. We 
consider the following class of vectors in the space T: 

$ = B {n}{m} ^ {n} ,{m} i 

(43) 

{n},{m} 

where -B{„}{ m } depend on the sets of numbers {n}, {m} so that 

B{ n }{m\ = 

if at least one of the conditions 

oo oo 

n a < aM, ]T m c < a/c ( 44 ) 

<J=1 c=l 

is not satisfied. The metastable state corresponding to our function system exists if the 
functional 

' <[>. ( H - OS ) <v 

- (45) 



has a minimum on the class of vectors (43) for some value of a. In the number of cluster 
types is > 2, the class of vectors is determined similar to (43) by using inequalities similar to 
(47). And the existence of a metastable state for a larger number of clusters is determined by 
the existence of a minimum of the corresponding functional of type (45) on the corresponding 
vector class. 
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We proposed a quantization of entropy and of free energy for many particles. 
The results can readily be carried over to the case of many fields. 

In particular, we can define quantum thermodynamics in the case of many fields by 
generalizing our operation of quantization over clusters to the case of third quantization [9] 
(this operation was performed for the second quantization). 

Let us perform the quantization of thermodynamics in the theory of many scalar fields. 
In [9], a third quantization of the model of N scalar fields is presented for the Hamiltonian 
of the model of the form 

— 1 N r N r 

H N = 2 J d0C {jt a (x)Tt a (x) + V(fa(x)V(fa(x) + VT? ' (p a (x)(f a (x)^ + E ^ J dx(pl(x)(f 2 b (x) , 
a=l a, 6=1 

(46) 

where x G R 3 , e, and m are given quantities and 7f a (a;) and (p a {.x) are operators in the space 
Cn of the symmetric functionals, ^jvfy?!^), . . . , <Pn(-)], of iV functions y?i(-)? • • • , Vn^-), 

(Pa(-)^N[(pl(-), • • • , <Pn(-)] = <Pa(-)^N[<fl(-), • • • , <Pn(-)], 

n a (-)VN[<Pi(-), <Pn(-)] = -i-r — tt^jv[^i(-)> • • • , <Pn(-)]- 
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Let us consider the Fock space C = ©^ =0 £jv whose elements are sequences of symmet- 
ric functionals \1/ A r[(/9 1 (-), . . . , ip N (-)], N = 0,1,..., and introduce operators of generalized 
functionals [</?(•)] on C as follows: 

(/ D V A+l^-)]X[ V (-)]^ ^ . . . , m-)] = 

1 k 

= -7= E^ba(0]**-lbl(0> ■ • ■ , Va-l(0, v«+i(0, • • • , ¥>*(0], ( 47 ) 

V fc a=l 

(/ AM>(.)]x>(-)]*) fe x MO, ■ ■ ■ , Vfc-i(-)] = 

= VkJ D<pX*[<p(-)]* k [<p(.), ^(0, ... , 

To the Hamiltonian of iV fields of the form (46), we assign the following operator H on the 
Fock space C: 

(#*)ivM-), • • • , MO] = #J\r*ivM0, • • • , MO], 
which can be represented as 

# = £ / £><M~>(0] / ^ (-^y^) + Vy?(:r)Vy2(:r) + mV(i)^)j ^b(01 + 

ityity I dV(^)0 2 (^)^ + b(-)]^"b(-)]^ + [0(-)]^"[0(-)]- (48) 

The expression (48) is the triply quantized representation of the Hamiltonian (46). 

Now we introduce the cluster spaces M.N u ...,N k , where k > 1, N±, . . . , N k > are integers. 
The elements of these spaces are functionals of the form 

*jvi,...,.^1(0>JiV" 

...;¥>*(•),..., ^(.), Ji, • • • A fc _ fc+ i(0, • • • > ^7v fe ("), JfcTvJ, (49) 
that are symmetric with respect to the permutations of 

V 9 pi+l(')) • • • > Vpi+i(')) Jp+15 V?gZ+l(0> • • • 5 V?gZ+j(0>.7g+l 

for all Z = 1, . . . , k and p, q = 0, . . . , iVj — 1, where j^, s = 1, . . . , JVj, take values in the set 
0, 1, . . . , We also introduce the cluster Fock space .M = ©5v n =o-M.N 1 ,...,N k , which con- 
sists of sequences of functionals (49). Similarly to (47), we can introduce operator generalized 
functionals Bf[<p[(-), . . .,$(■),?]. 

We introduce an operator measure p, which is an operator on C and depends on the 
functionals Bi[ip[(-), . . . , ip\(-), j 1 }: 

p = exp {-j D V h [<p(.)] A- [v^(-)]) x 

k h 

x exp (E E / • • • / MM . . . Dy\D4>\ 

l=lj'=0 J J 

xBt[^(.),...M(-)J]BtM-),...M(-)J]^ 

x i+ [^(.)] i- [^(oi i + M(-)] i- [<d(oi), (50) 
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where the numbers over the operators stand for the order of their execution [5]. We omit 
the arguments p in what follows. Furthermore, let 

H{Bt[.],B 1 [.],...B' i [-],B k [.]) = 

= ex P [EE /••• / D<p[... D<p\BM(.), . . . , rfCOlBMO), • • • , ¥>{(■)] ) x 
x Sp (pH) . (51) 



By the cluster representation of the triply quantized Hamiltonian fl48 ) we mean the 
following operator on A4: 

n = n (it [•], 4r [•],-■■ 4 [■], 4 (■]) • (52) 

We introduce the operator of entropy of many fields on Ai by the relation 

s = s reg (it [■}, ii [■},■■■ k [■}, k [■]) , (53) 

where S Teg is the regularized functional 

s = Sp ( ?ln fe))' (54) 

Definition. The numbers A such that 

(n + es)$ = \$, <$>eM, $^0, (55) 
are called quantum values of free energy of many fields under the temperature 6. 

In the general case of a system of many fields, the quantization of entropy and free energy 
can be carried out just as was done in the present note. 

9 Supersecondary quantization and quantization of en- 
tropy with preserved charge 

Earlier the supersecondary quantization and quantization of thermodynamics was introduced 
in the case in which the secondary quantized Hamiltonian of the system commutes with 
the operator of the number of particles. In this section the secondary quantization and 
quantization of thermodynamics with number of cluster types equal to 2 is generalized to 
the case of charge preservation. 

Let ni,n 2 ,m > be integers. Consider the functions 

^n,mi ,mi C^' 1 : Jlj • • • j *^ni i 3n\ i V\ j ^1 j • • • > Vni j ^ri2 1 ^1 1 ^1 > • • • ) ^"m > ^-'ro) (56) 
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that are symmetric with respect to permutations of the pairs x p ,j p and x q ,j q , with respect 
to permutations of the pairs y r ,i r and y s ,i s , and with respect to permutations of the pairs 
zi, wi and z t , w t , where the variables x p , y r , z t , wi take values on the three-dimensional torus 
T 3 , and the variables j p and i r attain values 1,2, .. . and are called numbers. The Hilbert 
space of functions (56) with the norm 

ll<3> II 2 = 

oo oo oo oo 

= ■■■ Yl ■ ■ ■ Yl ■ ■ ■ Y ■■■ dxt... dx ni dyi . . . dy n2 dz x dw\ . . . dz m dw m x 

il=l 3n]=y 11=1 in 2 =l 

X | ( J > n,mi,m,2 (-^l , Jl , • • • , ^ni i jni , 2/l, H , . . . , 2/n 2 , ^"2 > ^1 5 ^1) • • • j z mi | (57) 

we denote by T n ^n 2 ,m- Now we introduce the cluster space JF: 

oo 

•T 7 = ® rei,ra 2 ,mj (58) 

rii,n2,m=0 

whose elements are sets of functions (56): 

$ = {^ ni ,n 2 ,m}, n 1 ,n 2 ,m = 0, 1, . . . , 

oo 

||<£>|| 2 = "V II* II 2 
ll^ll ll^ni,n 2 ,m|| • 

ni,n2,m=0 

The space is a special case of the Fock space in which we introduce the creation and 
annihilation operators bf(x,j), bf(y,i), B ± (z, w): 

(Pi J ) ( &)ni,n2,m(%l , jl , ■ ■ ■ i %ni , Jni i 2/l, il , • • • , 2/n2 ? ^n 2 ? ^1 ? ^1? • • • j ^mi Wm) 

(iCl, Jl, • • • 



yi, ii, . . . ,t/ n2 ,« n2 ,zi,wi, . . . , z m ,w m J, 

(^1 J ) < ^ > )ni,ra 2 ,m(2'l ; Jl j • • • , X ni , J ni , J/i, %i , . . . , 2/ n2 , i n2 , Zi , "U7i, . . . , Z m , Wm) 



V"H^-^rii+l,ii2,m(''') J, -^l , Jl, ■ ■ ■ i X ni , J ni , 2/l, ili • • • i Vri2 1 ^n 2 1 z l 1 ^1, • • • 1 z mi ^m) 1 
(^2" (?/) *)^)iii"2,wi(*''l) Jl, • • • 5 *^rii 5 Jni , J/l, • • • , 2/n 2 , ^n 2 , ^1) • • • , w m) = 

X^l^ ~~ yi)®m,n2-l,m(xi,jl, ■ ■ ■,Xn 1 ,3n 1 ,yi,il, ■ ■ - , 
n 2 J=l 



y J+ i,i /+ i, . . . ,y n2 ,tn2, z i, w i, ■ ■ ■,z m ,w m ), 

(p2 (Vi i)^)ni,n2,m(Xl, jl, ■ ■ ■ , X ni ,j ni , J/i, ?i, . . . , |/n 2 ; ^n 2 ; ^1, Wlj • • • ; z mi U! m ) = 

\Al2 ~i~ l^Tii,Ti2+l,mO'l) • • • 5 -^ni 5 Jni 5 2/? ^, 2/l 5 • • • 5 2/n 2 5 ^n 2 ? z l 1 ^li ■ ■ ■ 1 z mi ^m) 5 

ni ,n2,mV*^l 5 Jl ; • • • ; X ni , J ni , yi , %i , . . . , ?/ n2 , %n2 1 z li^li • • • 1 z mi 
■y m 

,n 2 ,m—l (^lj Jl, • • • 5 X ni , jVii , J/l, i\i • • • , 2/n2 , ^n2 , 



Zi, Wi, . . . , Z^, Wi-x, Z l+1 ,W l+1 , ...,Z m , W m ), 
(B (z, 'U7)$) rilin2jm (Xi , ji , . . . , X ni , J n i 1 2/l, • • • , 2/«2 ) ^n2 5 ^1 5 ^1 5 • • • 5 Z mi W m ) 

= \Jv(l + l*ni,n2,m+l (>^1, Jl, • • • , -^ni , Jni , 2/l, ^1, • • • , 2/n 2 , ^n2 , ^ , ^ , ^1 , ^1, • • • 5 z mi W m ) , 
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where 8jf is the Kronecker delta. The operators bf(x,j) and 6 X (x,j) will be called the 

creation and annihilation operators for charge 1, the operators b^x.j) and will be 

- — - 

called the creation and annihilation operators for charge — 1, and the operators ' (x,j) and 
bi(x,j) will be called the creation and annihilation operators for charge 0. We consider 
the subspace of the space 0~ni,n2,mi which consists of functions ( |56|) symmetric with respect 
to permutations of the variables xi and x p , symmetric with respect to permutations of the 
variables x\ and z r , symmetric with respect to permutations of the variables y q and y s , and 
symmetric with respect to permutations of the variables y q and w t and denote this subspace 
by m . The following subspace of T: 



•F — ® F ni,n 2 ,m (59) 



ni ,ri2,m 

ni,n2,m=0 



will be called the boson subspace of the cluster space. 

Now we consider the "usual" Fock space 7i of boson systems with charges ±1, which 
consists of sets of the functions = {<f> qi ,q 2 (%i, yq 2 )}i <?i; <?2 = 0, 1, . . . sym- 

metric with respect to permutations of the variables x\ and x v and symmetric with respect 
to permutations of the variables y^ and y r . In this space the creation and annihilation oper- 
ators ip^(x) for bosons with charge 1 and the creation and annihilation operators ijjf(y) for 
bosons with charge —1 are defined as 

(tpf(x)(f)) quq2 (x u ...,x qi ,yi,..., y q2 ) = 
1 91 

22S(x - :r/)</> gi _i )g2 (zi, . . . x m , . . .,x qi ,y u . . .,y q2 ), 
<?i i=i 



($i(x)<f>) qi , q2 (xi, ■ ■ -,x qi ,yi, . . .,y q2 ) = 

= \jqi + 10 9l +i, 92 (a;,xi, . . .,x qi ,yt, . . .,y q2 ), 

(^d/)0)< ?1 ,g 2 ( a; i ) • • • ,x qi ,y u . . .,y q2 ) = 
1 91 

= z2s(y - yi)(Pqi,q2-i( x ii ■ ■ -i x qiiV^ ■ ■ -,yi-i,yi+i, ■ ■ -,yq 2 ), 

92 Z=l 



$2 {y)<f>)qi,q 2 (xi, ■■■,x qi ,yi,..., y q2 ) = 

= \Jq2 + 10 9ll92+ i(xi, . . .,x qi ,y,yt, . . .,y q2 ). 
An operator in Ti of the form 

Q = dxipi (x)ijji (x) — / dx ^ (x)^ ( x ) (60) 



is called the charge operator Q. Let A be an operator in the space 7i commuting with the 
charge operator Q. The supers econdary quantized operator A is the following operator in T: 

2 12 12 1 

A = A$t (x,j),bt (y,i),6j (y,i),B + (z,w),B~ (z,w)}, (61) 

where the numbers over the operators stand for order of their action and the functional 
A[6J(0, &i(0, &S(0, B*(-), B(-)} has the form 

A[bl(-), &£(.), b 2 (-), £*(•), 5(.)]= Sp (Ap[bt(-), MO, MO, MO, s *(0, £(0L) x 
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/ oo . oo . \ 

xexp -Ey dxblfafihixj)-^ J dyb* 2 (y,i)b 2 (y,i)\ , (62) 



x 



j=l i=l 

where &i(-)> ^(Oj B*(-), B(-)] is an operator valued functional ranging in the set 

of operators in the space Ti and having the form 

= exp Q J J dz x dw x B(zx, Wi) ^ {z x ) (wi)j x 

x exp (]^J J dz 2 dw 2 B*(z 2 , w 2 ) $f (22) (ufe)^ x 

(00 j2 _i 

dxidx 2 6i(ari, j)6J(a:2, 3) ■tpf (xi) $f (x 2 ) 

x ex P d 2/i rf 2/2 62(2/1, i)^(2/2, i) $t (l/i) ^2" (jte)^ x 

x exp ^- y dx (x) $f (x) - J dyfy (y) 4> 2 (y)j . (63) 

To each element $ of the boson subspace JF B of the space JF, we bijectively assign the set 
)3qi)^i) ■ ■ ■ ■> ^2)}) 9i> 92 — 0, 1, . . ., with values in the space of the functions of 
the numbers 

(0(jl) • • • 5 J<ji j • • • j ^q2))ni,n2 • • • 5 *^ni iVli ■ ■ ■ 1 Vri'i) = 

= fiqi—ni,q2—ri2@{ n l Ql)^qi,q2,ni—qi jit ■ ■ ■ ? x qu3qn Hit • • • j 2/<72> ^<?2 ) 

x qi+liUq2+li ■ ■ ■ 1 X n\)Un2) ) (64) 

where #(•) is the Heaviside theta function. Hence, to the secondary quantized operator A, 
there naturally corresponds an operator A in the space T B such that 

(74$)g li g 2 ?Tl (Xx , 3li ■ ■ ■ i Xq 1 , jq 1 , J/l, l±, ■ ■ ■ , 2/<j 2 , 2g 2 , X^+i , 2/<j 2 +l j • • • j %q 1 +rm Dq2+m) 

= (^0(j'l? • • • h, ■ ■ ■ , iq 2 ))qi+m,q2+m(Xl, ■ ■ ■ , ^gi+m, J/l, • • • , Uq2+m)- (65) 

We have the following theorem. 

Theorem. TTie projection of the supers econdary quantized operator A on the boson 
subspace is equal to A. 

We introduce an operator of supersecondary quantized entropy S for systems with pre- 
served charge 

— 2 J J 2 J 2 1 

S = S reg [bf (x,j),bi (x,j),b% (y,i),b 2 (y,i),B + (z,w),B~ (z,w)\, (66) 
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where the functional S reg [b*('), &i(-), ^(Oj^Oj B*(-), B(-)] is equal to the regularized func- 
tional 

S[bl(-)M-)M)M),B*('),B(-)]= 

m(-),bi(-),b* 2 (-),H-),B*(-),B(.)} \ ^ 



xexp -£/ rf y fe 2(y^) fe 2(y,«)j • ( 67 ) 
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